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Ultrasonic  nonspecular  reflectivity  near  longitudinal  critical  angle 

Tran  D.  K.  Ngoc  and  Walter  G.  Mayer 

Physics  Department,  Georgetown  University,  Washington,  D.C.  20057 

A  numerical  integration  method  is  developed  to  determine  the  intensity  profile  of  an  ultrasonic 
beam  reflected  from  a  liquid-solid  interface  near  the  longitudinal  critical  angle.  The  profiles  are 
calculated  for  different  combinations  of  frequencies  and  beam  widths  with  the  angle  of  incidence 
being  varied  about  the  longitudinal  critical  angle  for  a  water-Plexiglas  interface.  These 
calculations  demonstrate  the  existence  of  nonspecular  reflectivity  near  this  particular  critical 
angle  and  provide  a  quantitative  description  of  its  basic  features.  Theoretical  results  and 
experimental  measurements  are  compared. 

PACS  numbers:  43.20.Fn,  43.3S.Pt 


I.  INTRODUCTION 

An  analytical  formulation  to  treat  the  reflection  of  a 
bounded  ultrasonic  beam  at  a  liquid-solid  interface  was  first 
established  by  Schoch.'  That  theory  predicts  a  lateral  dis¬ 
placement  of  the  reflected  beam  for  all  total  reflection  an¬ 
gles.  This  simplified  model  fails  to  account  for  various  ex¬ 
perimentally  observed  phenomena  at  the  Rayleigh  critical 
angle.  Making  use  of  the  existence  of  a  pole-zero  pair  in  the 
reflection  coefficient  at  the  Rayleigh  critical  angle,  Bertoni 
and  Tamir2  describe  an  analytical  approximation  model  that 
js  able  to  explain  all  important  features  of  the  nonspecular 
reflectivity  at  the  Rayleigh  angle.  Their  approximation 
method  for  a  liquid-solid  interface  has  been  subsequently 
extended  to  a  study  of  reflectivity  from  a  solid  plate  im¬ 
mersed  in  a  liquid,3  the  diffraction  effect  observed  away  from 
the  interface,4  and  the  influence  of  attenuation3  in  the  solid 
medium.  Recently,  the  same  Rayleigh-angle  nonspecular  re¬ 
flectivity  problem  was  investigated  by  means  of  a  numerical 
integration  method  which  produces  results  consistent  with 
those  of  Bertoni  and  Tamir,2  but  also  allows  one  to  calculate 
changes  in  the  profiles  of  the  reflected  beam  for  angles  of 
incidence  near  the  Rayleigh  critical  angle.4 

While  nonspecular  reflectivity  near  the  Rayleigh  criti¬ 
cal  angle  has  been  studied  extensively,  little  is  known  about 
the  possible  nonspecularly  reflected  phenomena  for  bound¬ 
ed  beams  incident  at  or  near  the  longitudinal  and  shear  criti¬ 
cal  aflgles.  Experimental  evidence  indicating  the  existence  of 
nonspecular  reflectivity  similar  to  that  at  the  Rayleigh  angle 
was  reported  by  Neubauer7  for  Mallory  1000  at  the  shear 
critical  angle  and  for  Pyrex  glass  at  the  longitudinal  critical 
angle. 

This  paper  is  intended  to  determine  the  reflected  beam 
profile  of  a  bounded  sound  beam  incident  at  and  near  the 
longitudinal  critical  angle  of  a  liquid-solid  interface. 

II.  EFFECT  OF  ABSORPTION  ON  PLANE  WAVE 
REFLECTION  COEFFICIENT 

A  convenient  form  of  the  infinite  plane  wave  amplitude 
reflection  coefficient  A  for  a  bquid-tolid  interface  is  given1  by 

-  -  (*?  - 1*2  ?  +  -pkfa/K 

(ki,  —  lk\f  +  Ak\K,Kd  +pkfa/K  ’  } 

where 


p  =  liquid  density /solid  density, 
k  =  co/v, 
kd  =  eo/vd  , 
ks  =  w/vs  , 
kx  =  k  sin  6,  , 

K=  (k2-kl)l/1. 


FI 0. 1.  MifDitiKfe  tad  phase  of  the  nAectxn  cocAdeot.  The  tbaorptive 
cmcs  are  dsaoMd  by  (tasked  curves  sad  the  tonka  eases  by  hU  carves.  P 
denotes  a  wrtw-Pkodjto  interface  far  a  frequency  of  2  MIfa  sad  Stems 
a  wstir  ilsislm  mil  tnwrtkoe  for  IS  MHz. 


7948  J.  Appl.  Phys.  10(12),  Ososmbsr  1079 


0021-00 79/70/1270494)4001. 10 


®  1980  American  Insttuis  of  Physics  7949 


TABLE  I.  Sound  velocities  and  density  ratios. 


Materials 

(m/sec) 

V \ 

(m/sec) 

P 

Stainless  steel 

$840 

3130 

0.126 

Plexiglas 

2740 

1350 

0.84 

Water 

1480 

1 

where  co  is  the  angular  frequency,  0,  is  the  incident  angle,  u  is 
the  sound  velocity  in  the  liquid,  and  vd  and  v,  are  the  veloci¬ 
ties  of  the  longitudinal  and  shear  waves  in  the  solid,  respec¬ 
tively.  The  x  direction  is  along  the  interface,  and  the  inci¬ 
dence  is  from  the  liquid  half-space. 

It  is  known  that  the  behavior  of  R  (kx )  is  critically  de¬ 
pendent  on  the  attenuation  constants  of  the  media,  especial¬ 
ly  at  the  Rayleigh  critical  angle,  where  the  so-called  “least- 
reflection  frequency”  is  normally  observed  for  many  solids.’ 
Near  the  longitudinal  and  shear  critical  angles,  attenuation 
only  slightly  affects  the  magnitude  and  the  phase  of  R  (kx). 
However,  it  should  be  noted  that  inclusion  of  attenuation 
into  the  expression  for  R  (kx )  does  smooth  out  the  discontin¬ 
uous  behavior  of  the  first  derivatives  of  the  magnitude  and 
the  phase  of  R  (kx).  In  the  complex  plane  of  kx,  inclusion  of 
attenuation  results  in  the  property  that  the  branch  points 
±  k,  ±kd,  and  ±  k,  of  R  (kx )  and  their  corresponding 
branch  cuts  are  then  located  off  the  real  axis.  This  shift  of 
position  of  the  branch  points  makes  a  numerical  integration 
method  possible. 

The  effects  due  to  attenuation  in  materials  can  be  incor¬ 
porated  into  R  (kx)  by  introducing  the  complex  wave  num¬ 
bers  k,  kd,  and  k,  as  defined  below.  Denoting  a  as  attenu¬ 
ation  per  wavelength  for  the  sound  wave  in  the  liquid,  and  ad 
and  a,  for  the  longitudinal  and  shear  waves  in  the  solid,  the 
complex  wave  numbers  can  be  written  as 

k  —  (<o/v)  (1  +  ia/2ir) , 

kd  =  (o/vd)  (1  +  iad/2ir) ,  (2) 

k,  =  (ot/v.)  (1  -(-  ia,/2ir) . 

To  illustrate  the  influence  of  attenuation,  the  magni¬ 
tude  and  the  phase  of  R  (kx)  are  plotted  as  functions  of  the 
angle  of  incidence  for  water-stainless  steel  and  water-Plexig- 
las  interfaces  with  and  without  attenuation  (Fig.  1).  The  pa¬ 
rameter  values  used  to  produce  these  curves  and  all  subse¬ 
quent  figures  are  given  in  Tables  I  and  II. 

Figure  1  shows  that  the  curves  for  stainless  steel  are 
typical  of  many  other  solid/liquid  combinations  as  previous- 
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FIG.  2.  Reflected  beam  profiles  for  a  bounded  beam  of  0.75-in.  width  inci¬ 
dent  at  0L  onto  a  water-Plexiglas  interface.  The  dashed  line  represents  the 
incident  beam,  the  solid  lines  represent  the  reflected  beam  profile  for  (A)  1 
MHz,  (B)  2  MHz,  (C)  6  MHz. 

ly  reported,*  while  those  for  Plexiglas  exhibit  a  remarkably 
different  pattern,  expecially  near  the  longitudinal  critical  an¬ 
gle  (0L  =  32.7(F). 

III.  NONSPECULAR  REFLECTION  OF  A  BOUNDED 
BEAM  AT  LONGITUDINAL  CRITICAL  ANGLE 

Since  the  strongly  varying  behavior  of  R  (kx )  near  the 
Rayleigh  angle  leads  to  pronounced  nonspecular  reflectivity 
while  the  effect  is  insignificant  where  R  (kx  )  does  not  possess 
this  characteristic,*  one  may  expect,  on  the  basis  of  the  re¬ 
sults  shown  in  Fig.  I,  that  Plexiglas  may  also  show  nonspe¬ 
cular  reflection  for  incidence  near  the  longitudinal  critical 
angle. 

Using  the  method  of  representing  a  bounded  sound 
beam  by  a  Fourier  integral  in  k  space,2  a  numerical  integra¬ 
tion  algorithm  is  developed  on  the  basis  of  Simpson’s  Rule  to 
calculate  the  reflected  beam  profile  of  a  Gaussian  incident 
beam.  The  advantage  of  this  method  is  the  possibility  of  in¬ 
vestigating  nonspecular  reflectivity  for  any  angle  of  inci¬ 
dence  and  for  non-Gaussian  incident  beam  profiles  as  well  as 
the  Gaussian  one.  The  method  has  been  used*  to  study  the 
nonspecular  reflection  phenomena  near  the  Rayleigh  angle. 
Here  it  will  be  used  to  examine  reflectivity  in  the  angular 
range  about  the  longitudinal  critical  angle. 

It  was  shown  by  Bertoni  and  Tamir2  that  the  amplitude 
field  of  a  bounded  beam  reflected  from  a  liquid-solid  inter¬ 
face  can  be  described  by  an  integral,  which,  after  slight  modi¬ 
fications,  normalization,  and  several  changes  of  variables  to 
facilitate  the  integration  process,  is  given  by 

UrQC)  -  (2r)- f  +  W‘  R(Ax)V{Ax)EiAx*)dAx  ,  (3) 

where 

R(A  1.  (2Al-Api  +  ABdBrA l  - pA *Bd/B 
*x) "  (M  J  -  A  ))'  +  4 B4BJ  l  +pA',Bd/B ' 
V(4,)-Cap (4.  -A,?), 
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FIG.  3.  Reflected  beam  profiles  for  a  bounded  beam  of  frequency  2  MHz 
incident  at  9 ,  onto  a  water-Plexiglas  boundary.  The  beam  width  is  (A)  0.3 
in..  (B)  0.73  in.,  (C)  1.0  in.  The  dashed  line  represents  the  incident  beam. 


E(AxfX)  =»  txp(iCAxX) , 
where 

A  —  1  4-  ia/2ir , 

Ad  =  (v/vd)(.l  +  iad/2if) , 

A,  =  (o/u,Xl  +  ia,/lir) , 

B  =  (A1  —  A\)'n , 

Bd=(Al-Al)'/2, 

B, =(A2-Aiy'\ 

A,  =  sin  8,  , 

X  —  x/Wq  , 

C  =  6>wp/v , 

Vk=ir/C, 

with  x  being  the  coordinate  along  the  interface,  6,  the  inci¬ 
dent  angle,  and  2w0  the  beam  width  projected  onto  the 
interface. 

As  expected  from  observing  the  behavior  of  R  (kx)  for 
incidence  near  0L ,  the  calculated  reflected  beam  profiles  for 
the  case  of  water-stainless  steel  only  show  a  small  beam  dis¬ 
placement  of  about  2%  of  2jo0  which  is  normally  not  observ¬ 
able.  In  contrast,  for  the  case  of  water-Plexiglas,  the  dis¬ 
placement  of  the  intensity  peak  in  the  reflected  beam 
becomes  much  greater  and  can  be  observed  experimentally. 
This  characteristic  of  nonspecular  reflection  near  dL  is  not 
only  a  strong  (unction  of  the  material  properties,  but  also  of 
the  sound  frequency,  hence  of  ad  and  a„  the  beam  width, 
and  the  incident  angle. 

In  Fig.  2  a  set  of  beam  profiles  are  plotted  for  different 
frequencies  while  keeping  the  same  incident  beam  width. 
The  beam  displacement  is  found  to  be  larger  for  lower  fre¬ 
quencies.  Figure  3  showi  another  set  of  reflected  beam  pro¬ 
files  for  a  sound  frequency  of  2  MHz  with  the  beam  width 
0.S,  0.7$,  and  1.0  in. .  The  calculations  point  to  the  fact  that 
as  the  beam  width  is  decreased,  the  beam  displacement  be¬ 
comes  larger  but  the  peak  intensity  is  reduced. 
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(C) 


<D> 

FIG.  4.  Schlieren  photographs  showing  beam  displacement  and  changes  in 
sound  intensity  near  0L .  The  incident  angle  is  (A)  9,  -  T,  (B)  6,  -  I*.  (C) 
+  I*. 


IV.  VARIATION  OF  BEAM  PROFILE  AS  A  FUNCTION  OF 
ANGLE  OF  INCIDENCE  NEAR  dL 

Reflected  beam  profiles  change  as  a  function  of  the  an¬ 
gle  of  incidence.  Their  variation  can  be  characterized  by 
changes  in  the  beam  displacement  and  the  intensity  peak, 
which  are  only  substantia]  near  the  longitudinal  critical  an¬ 
gle  for  a  water-Plexiglas  interface.  The  numerical  integra¬ 
tion  method  is  now  used  to  determine  the  reflection  profiles 
of  a  beam  of  2-MHz  frequency  and  0.75-in.  width  when  the 
angle  of  incidence  is  varied  by  steps  of  0.5*.  The  same  phe¬ 
nomenon  is  also  monitored  by  using  the  Schlieren  visualiza¬ 
tion  technique. 

Figure  4  consists  of  a  series  of  Schlieren  photographs 
showing  the  gradual  change  in  profiles  of  the  beam  reflected 
at  angles  first  smaller  than  0L  and  becoming  larger  until 
exceeding  0L.  A  sharp  increase  in  the  beam  displacement 
and  the  reflected  intensity  can  be  detected  for  an  incident 
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FIG.  S.  Reflected  beam  proflJee  for  e  bounded  beam  of  frequency  2  MHz 
and  beam  width  of  0.75  in.  incident  at  a  water-Plexiglee  interim  at  (A) 
8t  +r,qB)0L+T,(Q»L  +  1*.  (D)  9l.  The  dashed  Mae  rrpreeeaa  the 
Incident  beam. 
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FIG.  6.  Peak  intensity  and  displacement  of  a  bounded  beam  between 
6t  —  2*  and  $L  +  }’  for  a  water-Plexiglas  interface,  beam  width  of  0. 75  in., 
and  frequency  of  2  MHz.  Curve  (A)  shows  the  theoretical  peak  intensity 
variation.  The  theoretical  and  experimentally  measured  beam  displace¬ 
ments  are  denoted  by  curves  (B)  and  (C),  respectively. 


angle  slightly  less  than  0L .  A  number  of  these  profiles  are 
theoretically  determined  and  shown  in  Fig.  5  for  dL  <8, 

+  3*. 

In  order  to  establish  the  exact  relationship  between  the 
beam  displacement  and  8, ,  calculations  are  done  for  smaller 
increments  of  the  incident  angle.  The  results  plotted  in  Fig.  6 
indicate  the  existence  of  such  a  sharp  increase  simultaneous¬ 
ly  in  both  beam  displacement  and  reflection  intensity.  Com¬ 
parison  of  Figs.  6  and  4  shows  qualitative  agreement  between 
theory  and  experiment.  Measurements  of  the  beam  displace¬ 
ment  from  the  Schheren  photographs  are  also  shown  in  Fig. 
6.  Both  calculated  and  experimental  curves  exhibit  the  same 
general  variation  pattern  for  the  beam  displacement  with 
respect  to  the  incident  angle.  The  uncertainty  in  the  absolute 
values  of  the  beam  displacement  is  due  to  the  fact  that  small 
differences  between  the  values  of  vd,  v„p,  ad,  and  a,  listed  in 


Tables  I  and  II  and  those  which  apply  to  a  particular  Plexi¬ 
glas  sample  will  cause  the  calculated  displacement  to  differ 
from  the  observed  one.  The  greatest  uncertainty  for  this  cal¬ 
culation  lies  in  the  values  of  the  shear-wave  velocity  and 
absorption  for  Plexiglas.  These  values  are  not  known  to  a 
high  degree  of  accuracy. 

V.  CONCLUSION 

The  above  calculated  results  obtained  by  a  numerical 
integration  method  and  experimental  measurements  estab¬ 
lish  the  fact  that  there  exists  nonspecular  reflectivity  near  the 
longitudinal  critical  angle  for  a  liquid-solid  interface.  This 
phenomenon  is  found  to  be  significant  only  for  certain  liq¬ 
uid-solid  combinations,  one  of  which  is  water-Plexiglas.  For 
incidence  near  the  longitudinal  critical  angle,  nonspecular 
reflectivity  is  most  prominently  characterized  by  a  substan¬ 
tial  beam  displacement  which  varies  with  the  incident  angle, 
the  beam  width,  and  the  frequency-dependent  attenuation  in 
materials.  In  general,  the  numerical  integration  method  is 
shown  to  be  applicable  at  a  range  of  incident  angles,  where 
no  successful  analytical  formulation  had  been  available. 
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Numerical  integration  method  for  reflected  beam  profiles 
near  Rayleigh  angle 

Tran  D.  K.  Ngoc  and  Walter  G.  Mayer 

Physics  Department,  Georgetown  University,  Washington,  D.C  20057 
(Received  16  July  1979;  accepted  for  publication  31  December  1979) 

A  numerical  integration  method  is  developed  to  calculate  the  intensity  profile  of  an  ultrasonic  beam 
reflected  from  a  liquid-solid  interface.  This  numerical  treatment  is  used  to  calculate  nonspecular 
reflectivity  at  a  range  of  angles  of  incidence  near  and  at  the  Rayleigh  angle.  Calculations  for  a 
water-stainless  steel  interface  are  compared  to  a  known  approximate  analysis  for  various  beam  widths 
and  frequencies.  The  theoretical  predictions  of  the  reflected  beam  profile  near  the  Rayleigh  angle  of 
incidence  are  compared  to  experimental  results. 

PACS  numbers:  43.20.Fn,  43.20.Rz 


INTRODUCTION 

Nonspecular  reflection  effects  for  ultrasonic  waves 
were  originally  predicted  and  experimentally  verified 
by  Schoch1  and  later  by  Neubauer2  and  Breazeale  et  al3 
for  bounded  beams  incident  on  a  liquid— solid  interface. 
The  simple  model  presented  by  Schoch  does  not  account 
completely  for  different  experimentally  observed  pheno¬ 
mena  at  Rayleigh  critical  angle,  such  as  the  point  of 
intensity  minimum  and  the  trailing  field.  These  dis¬ 
crepancies  have  been  eliminated  by  Bertoni  and  Tamir4 
who  use  a  complex  Laurent  series  expansion  to  simpli¬ 
fy  the  mathematically  complicated  expression  of  reflec¬ 
tion  coefficient  in  terms  of  its  complex  pole  and  zero. 
Bertoni  and  Tamir’s  approximation  gives  excellent  re¬ 
sults  at  the  Rayleigh  critical  angle  which  is  associated 
with  the  complex  pole  and  zero. 

Attenuation  in  the  media  involved,  when  included  in 
these  analyses,  further  complicates  the  physical  prob¬ 
lem.  Becker  and  Richardson1  pointed  out  the  different 
behavior  of  the  reflecticn  coefficient  with  absorption  be¬ 
ing  taken  into  account.  In  such  cases,  total  reflection 
is  not  expected  at  any  angle  of  incidence  and  the  reflec¬ 
tion  amplitude  and  phase  are  shown  to  deviate  signifi¬ 
cantly  from  the  lossless  case. 

Recently,  Bertoni  and  Hou*  have  modified  their  form¬ 
ulation  to  incorporate  attenuation  into  their  theoretical 
investigation  of  nonspecular  reflectivity.  Pitts  et  al.1 
have  extended  Bertoni  and  Tamir’s  results  to  a  more 
complicated  case  of  the  liquid-solid-liquid  plate 
structure.  Breazeale  et  al*  used  Bertoni  and  Tamir’s 
analysis  to  investigate  the  reflected  beam  profiles  ob¬ 
served  at  large  distances  from  the  liquid-solid  inter¬ 
face. 


results  when  one  varies  one  of  the  physical  parameters, 
namely,  the  beamwidth,  the  sound  frequency,  and  the 
incident  angle. 

I.  MATHEMATICAL  REPRESENTATION  OF  BOUNDED 
BEAMS 

In  Fig.  1,  a  Gaussian  sound  beam  with  a  beamwidth 
2w  is  shown  to  be  incident  at  an  angle  6t  from  a  liquid 
onto  a  solid  medium.  The  amplitude  field  of  a  sound 
beam  with  finite  lateral  dimensions  can  be  expressed  as 
a  sum  of  infinte  plane  waves  via  a  Fourier  transform 
pair.4  An  incident  bounded  beam  can  then  be  repre¬ 
sented  by 

Ulu.(x,z)=(2vVlf  K(*it)exp[i(*xx+*t4)Jdft.l  (1) 

and 

V(kJ)=f~UlK(x,0)exp(-ikxx)dx,  (2) 

«o 

where  the  time  dependence  exp  has  been  sup¬ 

pressed,  kx  is  the  x  component  of  the  wave  vector, 

VlMj  is  the  Fourier  transform  of  the  incident  beam 
amplitude  at  z  =  0,  and  k,  is  defined  by 

with  k  being  the  propagation  wavenumber  of  the  sound 
beam  in  the  liquid. 

The  width  of  the  area  projected  by  such  an  incident 
beam  onto  the  interface  is  2wa,  which  is  given  by 

w0=iosec8(. 

Thus  for  well-defined  beams,  the  range  of  integration 
in  Eqs.  (1)  and  (2)  effectively  is 


Bertoni  and  Tamir’s  approximation,  although  giving 
excellent  results,  has  some  limitations.  It  is  not 
readily  applicable  to  those  critical  angles  other  than 
the  Rayleigh  critical  angle  and  presents  considerable 
analytical  problems  when  extended  to  non-Gaussian 
beam  profiles. 

In  the  following,  the  reflected  profile  is  determined 
by  numerically  integrating  the  Fourier  integral  des¬ 
cribing  the  reflected  beam  amplitude,  taking  into  ac¬ 
count  the  attenuation  in  the  materials  involved.  The  re¬ 
flected  profiles  are  calculated  and  compared  to  known 


(*,  -  w/w0)  *  ht « (k,  +  x/w0) 


(3) 


FIG.  1.  A  bounded  sound 
beam  Incident  onto  the 
plane  interface  with  a  sol¬ 
id.  The  scale  of  the  x 
axis  is  chosen  arbitrarily 
for  the  center  of  the  in¬ 
cident  beam  to  be  at  100. 
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and 


(4) 


Briefly,  the  Bertoni  and  Tamir’s  formulation  requires 
that  the  pole  and  zero  locations  be  identified  before  the 
reflected  profile  can  be  determined  by  Eqs.  (8)  and  (9). 


-«e0 
where 

kt  =fesinfy. 

The  Fourier  integral  for  the  incident  beam  f4*  (x,  z) 
can  be  interpreted  as  being  a  superposition  of  an  infin¬ 
ite  number  of  incident  plane  waves  with  different  ampli¬ 
tudes  v(kM),  all  of  the  same  wavelength  but  incident  at 
different  angles  within  a  narrow  angular  range  about 
et.  Extending  this  interpretation,  one  can  represent 
the  reflected  beam  profile  by  Integrating  over  individual 
reflected  plane  waves  which  are  now  affected  by  the 
plane  wave  reflection  coefficient  R{kx)  as 

UR(x,z)=(2n)-1  f  R(kx) V(kI)exp[i(kxx  +*,*]<&,.  (5) 

At  the  interface  z  =0  this  expression  reduces  to 

UR(x,  0)=(2jr)'‘  j  rt(*„)K(*,)exp{tM)<*«.  (6) 

with  R(kx)  being  given  by® 


mk  ,  -  (*i-  2*|)2  +4k^KtK4  -Pk*jct/K 

{  *’ '  (k\  -  2 klf  +4k  fa  +f*X A  ' 


p  =  liquid  density /solid  density, 
k  =  w/v,  kt  =  u/vt ,  k,  =  u/v„ 

where  v  is  the  sound  velocity  in  the  liquid  and  v4  and 
v,  the  velocities  of  the  longitudinal  and  shear  waves  in 
the  solid,  respectively. 


II.  AN  APPROXIMATE  ANALYTICAL  SOLUTION  TO 
THE  REFLECTED  BEAM  INTEGRAL 

Expanding  R{kx)  in  a  Laurent  series,  Bertoni  and 
Tamir  were  able  to  derive  an  approximate  form  for 
R(k„)  which  yields  the  reflected  beam  profile  expressed 
by  Eq.  (6)  for  a  Gaussian  beam  incident  at  the  Rayleigh 
critical  angle.  With  the  incident  Gaussian  beam  being 
normalized  to  unity,  the  incident  and  reflected  beam 
intensity  profiles  can  be  written  as4 

0)1'  =exp[~2(x/woy]  (8) 

and 

0)| a  =  C4«(  *,  0>J » ,  (9) 

where 

R,  =(Re*>-  *0)/(Re  kp  -  kp) , 

=  1  -  [(x*71*#  to  *r/2)  exptr*)  erfc  (y)] , 

where 

ill,  =Pole  of  R(k  x) , 
kg  =  Zero  of  R(kt) , 
y =(w#Im  Vfc)  -  (x/te„) . 


III.  INCORPORATION  OF  ATTENUATION  IN  THE 
MEDIA 

Effects  of  attenuation  in  materials  can  be  incorporated 
in  the  wave  vectors  k,kt>  and  k,.  Denoting  a  ,  at,  and 
a,  as  attenuation  per  wavelength  for  the  sound  waves  in 
liquid,  the  longitudinal  waves,  and  the  shear  waves  in 
the  solid,  respectively,  one  can  write  the  wave  vectors 
for  the  absorptive  media  as 

k=(w/v)  (1  +ia/2v), 

*«  =  (w/r4)(l+M,/2»),  (10) 

k,  =  (uj/v,)(l+ia^2v). 

Inclusion  of  attenuation  moves  the  branch  points  k, 
kt,  and  k,  off  the  real  axis.  Analytically,  introduction 
of  attenuation  into  the  terms  contained  in  the  reflection 
coefficient  [Eq.  (7)]  makes  R(kt)  differentiable  even  at 
the  longitudinal  and  shear  critical  angles.  For  the  loss¬ 
less  case,  the  reflection  coefficient  is  frequency-inde¬ 
pendent.  In  practice,  this  is  not  the  case  since  one 
usually  observed  in  experiments  different  profiles  for 
different  frequencies. 

Bertoni  and  Hou®  show  that  the  values  of  kt  and  ka  will 
change  with  the  frequency-dependent  values  of  at  and 
a,,  making  the  simplifying  assumption  that  ad  -a,  and 
that  the  absorption  in  the  liquid  medium  can  be  ne¬ 
glected  for  a  water-stainless  steel  interface.  Under 
these  simplifying  conditions  the  value  of  k  in  Eq.  (10) 
reduces  to  a  real  quantity,  and  the  imaginary  parts 
are  the  same  for  both  kt  and  However,  if  one  uses 
distinct  values  for  at  and  as,  and  includes  attenuation 
in  the  liquid,  the  loci  of  k,  and  kQ  are  different  from 
those  obtained  under  these  simplifying  assumptions. 

The  results  of  such  calculations  for  a  water-stainless 
steel  interface  are  shown  in  Fig.  2(a).  The  values  of 
the  longitudinal  and  shear  wave  absorption  coefficients 
used  are  those  given  by  Fitch10  (at=  0.00017,  a, 

=  0.00024)  and  the  value  of  the  attenuation  in  water  was 


W 


FK3.  2.  Pole  and  zero  locations  at  R{kx)  calculated  for  dif¬ 
ferent  attenuation  parameters  in  the  frequency  range  1-40 
MHz  for  water-stainless  steel  Interface  when  (a)  «*•*«,  and 
«al)|b)((  -a,  and  a  «  0. 


tiso 


J.  Anoint.  Sot  Am.,  Vot  67,  No.  4,  April  1980 


T.  D.  K.  NgocendW.G.  Mayer:  Reflected  beam  profiles  UB0 


1.0 


assumed  to  be  a  =0.00026.  Figure  2(b),  which  is  taken 
in  part  from  Bertoni  and  Hou,®  is  reproduced  for  com¬ 
parison.  The  results  shown  in  Fig.  2(b)  are  based 
on  the  assumption  that  at  =a,  and  that  a  is  negligible. 
One  notes  that  the  loci  of  the  poles  and  zeros  in  Fig. 
2(a)  no  longer  lie  on  a  straight  line.  The  values  of 
Re(k,)/Re(k)  are  of  no  particular  importance  for  the 
present  discussion  since  they  depend  only  on  the  as¬ 
sumed  value  of  the  sound  velocity  in  water. 


IV.  BEAM  PROFILES  CALCULATED  BY  NUMERICAL 
INTEGRATION 

In  the  Bertoni  and  Tamir’s  derivation,  it  is  important 
to  note  that  one  must  identify  the  Rayleigh  pole  and 
zero  and  that  the  reflection  coefficient  is  approximated 
in  the  following  form: 

*(*,)=(*„-*„)/(*,-*,).  (11) 

It  is  shown®  that  Eq.  (11)  is  only  sufficiently  accurate 
when  the  incident  angle  is  at  the  Rayleigh  critical  angle. 
Away  from  the  Rayleigh  angle,  the  accuracy  of  Eq.  (11) 
is  reduced  rapidly.  It  is  realized  that  the  reflected  pro¬ 
file  expressed  by  Eq.  (6)  can  also  be  calculated  by 
numerical  integration.  Since  the  exact  form  of  the  re¬ 
flection  coefficient  is  to  be  used,  such  a  numerical  ap¬ 
proach  can  provide  an  accurate  profile  even  away  from 
the  Rayleigh  critical  angle  without  a  prior  knowledge 
of  the  pole-zero  location. 

After  normalization  with  respect  to  an  incident  beam 
with  a  Gaussian  profile  and  making  use  of  the  effective 
integral  limits  given  by  Eqs.  (3)  and  (4),  Eq.  (6)  can  be 
written  as 

Vr(  X)  =(2ff)_l/j  fAl'W*  R(A,)V(A,)E(X,Ax)  dAx,  (12) 
where 

RIA  1  <gd| -A\BjB 
(2A*  -  A2J4  +ABtB,A\  +a',b]/B  ’ 

V(A.)  =  Cexp{-(  C/2)’(A,  -  A,)2] , 

E(X,Ax)=exp(iCAxX), 


FREQUENCY 
A  1  MHz. 

B  2 MHz 
C  4MHz 


X(V^/50) 


FIG.  4.  Reflected  beam  profiles  for  beamwidth  equal  to  20 
mm  at  various  frequencies  using  both  Bertoni  and  Tamir's 
analysis  and  numerical  integration  with  for  a  water- 

stainless  steel  interface.  Solid  curves  are  those  calculated  by 
numerical  integration. 


follows: 

A  =1  +ia/2w,  At  =  (v/vt)(l+iat/2t), 

A,  =  (t’/ti)(l  4-  my 2  V) ,  B  =  (A2  -  A2)1'2 , 

Bt  =(AJ  -  A|)l/2 ,  B,  =  (A2 -Aj)l/2 , 

A(=sin0,,  X  =  x/w0,  C  -  umjJ  v  ,  W„=r/C. 

The  reflected  profiles  which  are  shown  in  the  next 
section  are  numerically  determined  according  to  Eq. 
(12). 


V.  CALCULATED  AND  EXPERIMENTAL  RESULTS 

In  this  section,  the  reflected  profile  is  calculated  by 
two  approaches  for  different  sets  of  physical  parame¬ 
ters.  First,  the  frequency  is  kept  constant  while  the 
beamwidth  is  being  changed.  Then,  the  beamwidth  is 
held  constant  as  one  varies  the  sound  frequency.  In 
these  two  cases,  the  beam  is  incident  at  the  Rayleigh 
angle  and  the  values  of  k,  and  k0  substituted  in  the  re¬ 
sults  derived  by  Bertoni  and  Tamir  are  obtained  from 
the  attenuation  data  which  are  entered  into  the  numeri¬ 
cal  integration  approach.  This  was  done  to  eliminate 
the  possible  discrepancy  which  may  arise  due  to  differ¬ 
ent  data  used  for  attenuation.  Finally,  the  profiles  of 
the  reflected  beam  is  calculated  when  the  incident  angle 
is  varied  near  the  Rayleigh  critical  angle. 


where  the  new  dimensionless  variables  are  defined  as  In  Fig.  3,  the  reflected  beam  profile  is  plotted  for 

different  beamwidths  at  2  MHz  for  a  water-stainless 


100  200  300 

X  (Vi(/50) 

FIG.  3.  Reflected  beam  profiles  at  2  MHz  for  various  beam- 
widths  using  both  Bertoni  and  Tamir’s  analysis  and  numerical 
integration  with  ©  ,»©r  for  a  water- stainless  steel  interface. 
Solid  curves  are  those  calculated  by  numerical  method. 


100  200  300 

XCWb/50) 

FIG.  5.  Reflected  beam  profiles  for  beamwidth  equal  to  20  mm 
at  2  MHz  using  numerical  Integration  when  6,  Is  slightly  var¬ 
ied  from  ©f  for  a  water-stainless  steel  Interface. 
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FIG.  6.  Schiieren  pictures  of  beam  profiles  at  five  different 
angles  of  incidence  about  0r  for  a  water-stainless  steel  inter¬ 
face.  A:  e(  =  er-l°;B:  G,  =  6,  -  0.5";  C:  6(  =  Gr;D:  6,  =  6r  +  0.5°; 
E:  G,  =  Gr  +  1.0°. 

steel  interface.  Results  obtained  by  both  methods, 
Bertoni  and  Tamir’s  analysis  and  the  numerical  inte¬ 
gration,  are  presented  in  this  figure.  Also  for  the 
water-stainless  steel  interface,  the  reflected  profiles 
are  plotted  in  Fig.  4  for  various  frequencies  while  the 
beamwidth  remains  unchanged.  Again,  both  results  are 
shown  in  the  same  figure.  It  is  seen  that  Bertoni  and 
Tamir’s  analysis,  although  using  an  approximate  form 
for  R(kx),  agrees  very  well  with  an  almost  exact  calcu¬ 
lation. 

Figure  5  shows  the  results  for  a  2-MHz  beam  inci¬ 
dent  on  a  water-stainless  steel  interface  when  the  in¬ 
cident  angle  is  varied  from  the  Rayleigh  critical  angle. 
A  gradual  change  in  the  beam  profile  is  noticeable,  in 
particular,  the  slow  increase  in  intensity  in  the  “null 
intensity”  as  the  specularly  and  nonspecularly  reflected 
sections  of  the  beam  merge  into  one  with  lessening 
beam  displacement. 

The  same  gradual  change  is  observed  experimentally. 
Figure  6  consists  of  a  series  of  Schiieren  photographs 
of  the  incident  and  reflected  beams  for  a  water-stain¬ 
less  steel  interface.  The  beamwidth  and  the  sound  fre¬ 
quency  used  we  re.  19  mm  and  1.96  MHz.  The  incident 
angle  was  varied  by  steps  of  0.25°. 

It  should  be  noted  that  the  profile  of  the  ultrasonic 
beam  used  was  nearly  Gaussian  so  that  agreement  be¬ 
tween  theory  and  experiment  exists  qualitatively.  It 
should  also  be  noted  that  the  same  type  of  agreement 
exists  between  experiment  and  the  predictions  of  the 
present  theory  for  reflection  at  angles  greater  or 


smaller  than  the  Rayleigh  angle  by  up  to  2°.  This  range 
is  sufficiently  large  so  that  differences  between  the  ap¬ 
proximate  formulation  and  the  numerical  integration 
method  become  noticeable. 

VI.  CONCLUSION 

It  is  shown  that  a  numerical  integration  method  can 
be  used  to  determine  accurately  the  reflected  beam  pro¬ 
file  for  a  liquid-solid  interface.  The  method  is  not 
restricted  to  calculations  at  the  Rayleigh  angle  but  can 
be  extended  to  other  angles  of  incidence,11  where  the 
analytical  analysis4  cannot  be  applied.  The  results  ob¬ 
tained  by  the  two  methods  do  agree  at  and  very  near 
the  Rayleigh  critical  angle.  Schiieren  photographs  of 
the  reflected  beam  at  and  near  the  Rayleigh  critical 
angle  show  qualitative  agreement  between  the  results  of 
the  numerical  integration  method  and  experiments. 
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NONSPECULAR  TRANSMISSION  EFFECTS  FOR  ULTRASONIC  BEAMS 


INCIDENT  ON  A  SOLID  PLATE  IN  A  LIQUID 


K.W.  Ng,  T.D.K.  Ngoc 
J.A.  McClure,  W.G.  Mayer 
Department  of  Physics 
Georgetown  University 
Washington,  DC  20057 


The  nonspecular  phenomena  are  investigated  theoretically  for  an  ultra¬ 
sonic  beam  transmitted  through  a  solid  plate  immersed  in  a  liquid.  The  ana¬ 
lytical  method  used  by  Bertoni  and  Tami r  [  3  1  to  describe  the  nonspecular 
profile  of  a  beam  reflected  from  a  liquid-solid  interface  is  extended  to  the 
problem  under  consideration,  taking  into  account  all  existing  poles.  To  solve 
the  integral  represeating  the  transmitted  beam,  the  ampl i tude -plane  wave  trans¬ 
mission  coefficient  is  replaced  by  a  simpler  approximate  form.  The  transmitted 
beam  profile  is  then  calculated  from  both  single-pole  and  multiple-pole  formu¬ 
lations. 


t.  Introduction 


Since  the  work  of  Goos  and  Hanchen  in  19^7,  it  has  been  known  that  light 
waves  are  not  always  reflected  geometrically  but  may  be  shifted  in  the  specular 
plane  at  some  angles  of  inc i dence  [  1 3 •  Similar  displacements  have  been  observed 
with  sound  beams  as  well  as  other  non-specular  phenomena,  namely,  a  null  strip 
in  the  reflected  beam,  and  a  trailing  sound  field  which  follows  the  primary 
reflection  and  decreases  in  intensity  as  it  moves  down  the  interface.  In  the 
early  investigations  of  these  effects,  Schoch  [2]  proposed  a  simple  theory  which 
explained  the  displacement  of  the  sound  beam  but  failed  to  account  for  other 
non-specular  features. 

More  recently  Bertoni  and  Tami r  [3]  have  proposed  a  theory  which  success¬ 
fully  accounts  for  all  observed  non-specular  phenomena  in  reflection  from  a 
liquid-solid  (L/'S)  interface.  These  authors  emphasize  the  use  of  complex  poles 
of  the  plane  wave  reflection  coefficient  in  determining  both  an  approximate 
mathematical  form  of  the  coefficient  itself  and  in  finding  an  accurate  descrip¬ 
tion  of  the  reflected  sound  intensity  profile.  Pitts  [4]  has  extended  the  approach 
of  Bertoni  and  Tami r  to  the  study  of  the  reflection  of  a  sound  beam  from  a 
1 iqui d-sol i d- I i qu i d  (l/S/L)  systen. 

Experimental  observation  has  shown  that  non-specular  phenomena  occur  in 
the  transmitted  beam  as  well.  In  this  investigation  we  will  take  the  general 
theoretical  approach  used  by  Bertoni  and  Tami r  and  by  Pitts  and  extend  it  to 
the  problem  of  the  transmission  of  a  bounded  acoustic  beam  through  an  L/S/L 
system  [5]. 

1 1 .  Theory  of  the  Reflection  and  Transmission  of  a  Bounded  Beam  by  a  Solid  Plate 
Consider  a  bounded  acoustic  beam  of  width  2W  incident  at  angle  0.  on  a 
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solid  plate  of  width  d  immersed  in  a  liquid.  The  coordinate  system  is  as  shown 
in  Fig.  1  where  the  incident,  reflected  and  transmitted  beams  are  depicted. 

The  plate's  incident  surface  is  taken  to  form  the  xy  plane.  The  sound  field 
is  assumed  to  be  uniform  in  the  direction  perpendicular  to  the  plane  of  incidence 
and,  therefore,  has  no  spatial  y  dependence.  Using  Fourier  integrals,  the 
particle  displacement  amplitude  associated  with  the  incident  sound 

field  can  be  expressed  as  a  sum  of  infinite  plane  waves. 
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Here  k  is  the  projection  of  the  incident  wave  vector  on  the  x-axis,  k  =  k  sin0, 
x  * 

2  2  2 

with  6  the  angle  of  incidence  of  the  plane  wave,  and  k^  =  k  -  k^  .  The 
amplitude  V(k^)  of  the  plane  wave  components  of  the  incident  beam  can  be  ob¬ 
tained  by  Fourier  inversion  from  the  initial  beam  on  the  z  =  0  or  incident 
surface  of  the  plate. 


°° 
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The  Fourier  integral  of  the  incident  bean  can  be  interpreted  as  being 
composed  of  an  infinite  number  of  plane  waves,  all  of  the  same  wavelength  but 
at  different  incident  angles  8  with  respect  to  the  normal  to  the  surface.  If 
this  interpretation  is  extended  to  the  sound  fields  reflected  and  transmitted 
by  the  plate,  they  too  can  be  expressed  as  a  sum  over  individual  plane  waves 
by  using  the  reflection  and  transmission  coefficients,  R(k  )  and  T(k  )  for 
the  L/S/L  system, 

oo 

ea 
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Any  evaluation  of  the  reflected  and  transmitted  sound  beams  depends,  among 
other  things,  on  the  amplitude  distribution  of  the  incident  beam  V  (H,  o') 

Most  transducers  produce  roughly  a  Gaussian  bean  for  which  a  reasonable  approxi¬ 


mation  is 


<c 
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where  W  is  the  half  beam  width  and  £  is  a  coordinate  transverse  to  the  beam 
propagation  direction.  A  convenient  form  for  the  incident  profile  is 


U(x,<0  --  ['  (^)  4  L  1  , 


where  W  A^C  ^  Qj 


,  and  6.  is  the  beam's  incident 


angle  as  defined  in  Fig.  1.  Combining  (4)  and  (2)  we  get 

V(&x)--  Wo  ]  (5 

and  for  the  reflected  and  transmitted  sound  field  amplitudes  at  the  appropriate 


surfaces 
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An  evaluation  of  these  integrals  and  their  comparison  with  experiment  is 
a  central  problem  in  the  study  of  the  L/S/L  system.  The  reflected  beam  has 
been  thoroughly  investigated  by  Pitts  and  others,  and  in  this  study  we  shall 
extend  their  treatment  to  the  transmitted  sound  field. 


The  Amplitude  Transmission  Coefficient 


The  reflection  and  transmission  coefficients  which  appear  in  (6)  and  (7) 


are  obtained  from  the  boundary  conditions  at  the  surfaces  of  the  plate.  A 


straightforward  but  laborious  calculation  yields 
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where  the  following  variables  have  been  defined: 

P  P 1  iqu  i  d^sol  i  d  , 

k  =  kcos0  =  k  ( 1  - s  i  n  6)  , 

Ks  =  kscoses  =  kI(v/vs)2  -  sin2e]^, 

2  2  ^ 

Kj  =  kjCose^  =  k [ (v/v^)  -  sin  0]  , 

P  =  dKd’ 

Q  =  die  . 


The  coefficients  (8)  and  (9)  satisfy  the  relationship 


•' 
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appropriate  to  a  lossless  medium. 

Because  of  the  complicated  dependence  of  R(k^)  and  T(k^)  upon  k^,  a  direct 
evaluation  of  the  profile  integrals  (6)  and  (7)  is  very  difficult.  Complex 
contour  integration  is  a  possibility  but  the  presence  of  branch  cuts  in  R  and 
T  poses  difficulties  if  this  approach  is  taken.  Bertoni  and  Tami r  in  the  case 
of  reflection  from  the  L/S  interface,  and  later  Pitts  treating  the  L/S/L  case, 
have  addressed  this  problem  by  using  an  approximation  for  R(k^).  The  approxi¬ 
mate  reflection  coefficient  of  the  general  form 


(14) 


is  constructed  to  have  the  same  poles  p. ,  p.  =  B.  +  in.,  and  the  same  zeros  z. 

J  J  J  J  J 

as  R(k  )  but  no  other  singularities.  With  R(k  )  replaced  by  R'(k  ),  the  integral 
(6)  can  be  evaluated  using  the  residue  theorem. 

We  will  follow  that  approach  in  this  investigation  and  will  construct 
approximations  to  the  transmission  coefficient.  With  these  approximations 
the  profile  of  the  transmitted  beam  (7)  will  be  evaluated  and  compared  with 
exper i ment . 

In  a  lossless  medium  the  sum  of  reflected  and  incident  energies  must  equal 
the  incident  energy,  a  condition  expressed  mathematically  in  Eq.  (13)-  In 
the  present  investigation,  where  we  shall  try  to  use  as  much  as  possible  the 
earlier  results  of  Pitts,  this  equation  will  be  the  basis  of  our  approximation 
procedure  for  T.  Equation  (13)  states,  among  other  things,  that  whenever |R I*  0, 

| T J  *  I.  Two  very  simple  forms  for  T  which  have  this  property  are 


05) 


06) 


and 
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in  which  z.  and  p.  are  respectively  the  zeroes  and  poles  of  R(k  )  already 
J  J  ■* 

introduced  in  connection  with  (14). 


Each  of  these  forms  has  its  shortcomings.  Let  us  denote  the  relative 
phase  of  T,  T*  ,  and  T‘ 1  by  ij/,  ip  *  ,  and  i(i‘  1  ,  respectively.  Then  numerical  eval¬ 
uation  of  (35),  (lb),  (9)  and  comparison  show  good  agreement  between  ( T 1 (  and 
| T |  but  poor  agreement  between  the  phases  \J>'  and  ij>.  On  the  other  hand,  there 
is  good  agreement  between  '  '  and  ijr,  but  poor  agreement  between  j T '  '  J  and  |T  |. 


The  Transmitted  Beam  Profile 

In  addition  to  the  complexities  of  T,  the  integral  (7)  for  the  transmitted 
beam  is  complicated  by  the  presence  of  the  factor  exp  (dkz)  in  the  integrand. 

This  factor  does  not  appear  in  (6)  because  of  the  choice  of  the  incident  surface 
of  the  plate  as  the  z  =  0  surface.  This  problem  can  be  addressed,  however, 
by  expanding  kz  in  a  Taylor  series  about  sinB. 

A**- ^  /j-fcoap  O  +•  • .  .(17) 

All  the  terms  are  now  in  the  proper  form  for  the  evaluation  of  the  trans¬ 
mitted  beam  profile  integral.  Because  of  the  complexity  of  the  expressions, 
it  is  worthwhile  defining  an  integral  which  plays  an  important  role  in  the 
evaluation.  It  describes  the  profile  of  the  incident  beam  at  the  z  =  d  plane 
as  if  there  were  no  solid  plate  present.  Substitution  of  Eqs.  (2)  and  (17)  into 
(1)  yields 

oo 
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After  some  manipulation  this  can  be  shown  to  become 


fcN  W 

U(V1  -- 

in  which  the  variables 


i-m 


4  t  t. 


x^.  =  xcosQ  -  zsinQ, 

=  xsin6  +  zcos6,  (2C 

wT  =  W{1  -  (2iz/kW2cos6)3!i, 

have  been  introduced.  Equations  (19)  and  (20)  not  only  describe  the  Gaussian 
incident  beam  in  terms  of  a  new  set  of  coordinates,  2^  parallel  with  and  x^.  trans¬ 
verse  to  the  direction  of  incidence,  but  also  explicitly  express  the  beam 
spreading  caused  by  the  geometrical  diffraction  of  a  bounded  beam. 

The  transmitted  sound  field  amplitude  (7)  may  now  be  evaluated  for  either 
of  the  approximations  (15)  or  (16).  The  procedure  is  essentially  the  same  for 
both,  and  we  will  present  it  here  only  for  T"(k  ).  First  noting  that 


T".-rf  fi.  (- 


■)  =  !  \  ''•'*-*3  \  y 

then  substituting  into  (7),  we  find 
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Tr‘°'  \  /  _ 
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Here  1^  is  the  integral 
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where 
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I'  '  can  be  evaluated  analytically  by  using  the  convolution  theorem  in  a 
manner  like  that  employed  by  Pitts  in  treating  the  reflected  beam.  We  first 
wr  i  te 


-o o 


in  which 


o£> 
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G  (b-y)  is  integrable  by  complex  contour  integration.  For  b  >  y,  the  contour 

must  enclose  those  poles  with  cij,...,a  >  0  and  for  the  b<  y  case,  it  must 

enclose  those  with  a.,...,a  <  0.  In  forming  both  the  approximations  T*  (k  ) 

i  N  x 

and  T,,(k  )  only  poles  in  the  first  quadrant  were  used.  Therefore 

/i\  t 

and  the  convolution  integral  can  be  carried  through  to  yield 

l  (£J  ^  > 
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where 


v  - 


#rp  ; 


J  a  o 


Finally  the  profile  of  the  transmitted  beam  for  the  approximation  T'1,  using 
N  poles,  is 

t,  [c^\t 

\*n  7  '  ' 

*  -67f  (0  -^5°  CO]. 


The  other  approximation  rp  can  wr't 


ten  as 
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To  express  t'^  in  terms  of  T*^,  Eq.  ( 1 6 )  is  rearranged  to  give 


Denot i ng  0-.  ~ 


■9. 


T"%,V- 


,  it  is  noted  from  Eqs.  (30)  and  (31)  that 
A/  VW) 

-L  T  (*,) 


Where  t'!  (k  )  is  defined  as  the  sum  of  all  possible  products  a  ...  a  ,  in 
N  x  nl 

which  n^s  can  assume  any  value  from  1  to  N  and  must  satisfy  the  requirement 

a  <  a  <  ...  <  a  .  For  example,  if  N  =  3,  then  T*'^  =  a  +  a_  +  a  , 
nl  n2  nM  3  3 

T-3(2)  =  ?ja2  +  aja3  +  a2a3’  anc*  ^3^  =  a|a2a3‘  Insertion  of  Eq.  (32)  into 
Eq.  (7)  would  lead  to  the  transmitted  beam  profile  corresponding  to  T  ^  (k  )  , 


IV  (XyaV-2L  V  <l(^ 

T  M  *  1  ij 


where  XJ  u(m)  can  be  derived  from  Eq.  (29). 

/v 


V.  Results  and  Conclusions 


Theoretical  predictions  for  these  angles  of  incidence  were  obtained  by 
numerically  evaluating  the  transmitted  beam  profile  U.^  as  given  by  Eq.  (29). 
The  approximation  T",  Eq.  (16)  depends  upon  the  complex  poles 
the  transmission  coefficient  and,  as  shown  by  Pitts  et  al,  k^  *  Bj  corresponds 
to  the  condition  for  the  existence  of  a  vibrational  mode  of  the  plate.  Since 
k  *  ksinO,  p.  =  Rsin8.  implies  that  at  an  angle  of  incidence  0.  the  mode  of 
vibration  corresponding  to  the  jth  pole  is  excited.  At  an  angle  of  incidence 
very  near  the  jth  pole,  i.e., 
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<  O.  0  3  , 

the  contribution  for  the  jth  pole  to  the  total  intensity  profile  is  five 
or  more  orders  of  magnitude  greater  than  the  contribution  from  other  poles. 

At  such  angles  the  jth  pole  term  dominates  the  entire  expression  for  the  trans¬ 
mission  profile. 

In  our  calculations,  (x,z)  is  used  since  poor  agreement  between  } T J 

and  |T"|  is  not  tolerable  in  the  determination  of  the  transmitted  profile. 

The  transmitted  profile  was  first  calculated  for  the  case  where  6.  =  6 . .  Be¬ 
cause  of  the  dominance  of  the  jth  pole,  only  that  pole  was  accounted  for  in 
u|^(x,z)  and  the  resulting  profiles  are  shown  in  Fig.  2.  The  input  values 
are  chosen  for  a  brass  plate, in  water  with  fd  =  2  MHz. mm  and  sinBj  =  0.3&  for 
two  frequencies,  1  MHz  and  7  MHz.  The  outstanding  feature  seen  here  is  a 
lateral  displacement  of  the  transmitted  beam  relative  to  the  incident  beam. 

It  is  noted  that  in  Fig.  2  the  7  MHz  profile  has  a  larger  beam  displacement 
but  with  a  much  lower  intensity  peak. 

When  there  existmony  poles  in  the  neighborhood  of  sine.,  all  the  involved 

(N) 

poles  must  be  included  in  the  evaluation  of  U^,  (x,z)  to  obtain  the  appropriate 
profiles.  For  example,  the  transmitted  profile  was  again  calculated  for  a 
water-brass-water  system  with  fd  =  6.5  MHz. mm,  f  =  2MHz,  and  sin0.  varied  from 
0.31  to  O.36.  In  this  angular  range,  there  are  two  poles.  The  results  of  these 
calculations  are  shown  graphically  in  Fig.  3*  The  profiles  exhibit  some  inter¬ 
ference  effect  and  are  seen  to  change  with  the  incident  angle. 

The  method  employed  in  this  work,  while  giving  some  reasonable  results 
is  limited  in  several  ways.  One  is  that  the  restraint  (13)  between  R  and  T 
is  bothersome  and  neither  T*  or  T",  the  forms  of  which  are  suggested  partly 
by  it,  are  fully  satisfactory.  A  more  serious  limitation  is  the  dependence 
of  this  method  on  the  location  of  the  poles  of  T;  the  closer  sinB.  is  to  a  part 
cular  Pj  the  better  it  works.  A  technique  which  does  not  involve  the  poles 
might  be  less  restricted. 
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Restrictions  on  the  Existence  of  Leaky  Rayleigh  Waves 

NEAL  G.  BROWER,  DOUGLAS  E.  H1MBERGER, 
and  WALTER  G.  MAYER 

Abstmct-The  Rayleigh  wave,  an  inhomogeneous  surface  wave,  exists 
for  all  isotropic  elastic  solid  infinite  half-spaces.  When  the  free  surface 
of  the  solid  is  bounded  by  a  liquid  a  leaky  Rayleigh  wave  does  not  nec¬ 
essarily  exist  for  all  liquid/isotropic  solid  systems.  The  well-known  con¬ 
dition  for  the  existence  of  a  leaky  Rayleigh  wave,  the  sound  velocity  in 
the  liquid  must  be  less  than  the  shear  wave  velocity  in  the  solid,  is 
shown  to  be  a  necessary  but  not  a  sufficient  condition.  Additional  con¬ 
ditions  on  density  ratios  and  velocity  ratios  are  given.  Examples  are 
listed  showing  liquht/sctid  combinations  which  satisfy  the  liquid-shear 
wave  velocity  condition  but  not  the  additional  restrictions  and  thus  do 
not  support  a  leaky  Rayleigh  wave. 


I.  Introduction 

The  Rayleigh  wave  is  a  plane  inhomogeneous  surface  wave 
which  can  propagate  undiminished  along  the  surface  of  an  elas¬ 
tic  half-space.  The  wave  amplitude  decreases  exponentially 
normal  to  the  surface.  The  existence  of  such  surface  waves, 
which  can  be  supported  by  all  elastic  media,  was  shown  by 
Lord  Rayleigh  11]. 

There  also  exists  surface  waves  for  the  case  of  a  liquid/elastic 
solid  system.  One  surface  wave,  the  leaky  Rayleigh  wave,  is  an 
inhomogeneous  wave  which  is  damped  along  the  interface 
formed  by  the  adjoining  two  media.  The  damping  is  a  mani¬ 
festation  of  radiation  into  the  liquid  along  the  path  of  propa¬ 
gation.  In  the  limit  where  the  density  of  the  liquid  tends  to 
zero,  the  leaky  Rayleigh  wave  becomes  a  free  Rayleigh  wave. 
Unlike  the  free  Rayleigh  wave,  the  leaky  Rayleigh  wave  does 
not  exist  for  all  real  combinations  of  liquid/elastic  solid  media. 
A  well-known  example  is  an  elastic  solid  in  which  the  bulk 
shear  wave  velocity  is  less  than  the  sonic  velocity  in  the  ad¬ 
joining  liquid. 
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The  purpose  of  this  paper  is  an  investigation  of  further  con¬ 
ditions  on  the  existence  of  a  leaky  Rayleigh  wave.  The  condi¬ 
tions  involve  restrictions  on  the  following  ratios  which  describe 
the  liquid/elastic  solid  system:  P/./Ps,  VtIVD.  V LIVS,  where 
pL  and  ps  are  the  densities  of  the  liquid  and  the  solid,  respec¬ 
tively,  and  VL,  VD,  and  V$  are  the  sound  velocity  in  the  liq¬ 
uid,  the  longitudinal  and  shear  bulk  wave  velocities  in  the 
solid. 

Limits  on  these  ratios  are  derived  from  numerical  solutions 
of  the  secular  equation  for  leaky  Rayleigh  waves  and  specific 
examples  are  given  of  liquid/elastic  solid  combinations  which 
exceed  the  limits. 

II.  Theoretical  Background 

Consider  an  unbounded  isotropic  bulk  medium.  In  this  case 
two  and  only  two  different  waves  can  be  propagated.  If  the 
medium  is  bounded,  other  waves  may  exist. 

A.  The  Rayleigh  Wave 

For  an  isotropic  elastic  solid  with  a  plane  boundary,  a  sur¬ 
face  wave  always  exists.  The  wave  is  known  as  a  Rayleigh 
wave.  Briefly,  one  seeks  a  solution  to  the  equations  of  motion 
for  an  elastic  solid  where  the  displacement  amplitude  decays 
exponentially  perpendicular  to  the  surface.  At  the  boundary, 
all  stress  components  must  vanish.  This  leads  to  a  secular 
equation  for  the  velocity  of  the  surface  wave  V.  The  secular 
equation  is  given  by 

4  {VSIV)2[\  -  (VsIVfVl'UVsIVo)2  -  {VSIV)2)X'2 

+  11  -  2(VS/V)2]2  =  0  (I) 

where  VD  and  are  the  longitudinal  and  shear  bulk  wave 
velocities,  respectively.  There  are  three  roo(ts  to  the  secular 
equation  (1).  The  range  of  the  Poisson  ratib  o  determines  the 
nature  of  the  roots.  The  Poisson  ratio  may  be  written  in  terms 
of  the  bulk  velocities  as 

o*(Vh-2V$)IUVh-  V$).  (2) 

If  the  Poisson  ratio  is  in  the  range  where 

a  <  0.263  •  •  •  (3) 

then  there  exist  three  real  roots  to  ( 1 ).  However,  if  the  Poisson 
ratio  falls  in  the  range  where 

o>  0.263  ■  •  •  (4) 


001 8-9537/79/0700*0306 $00.75  ©  1979  IEEE 


IEEE  f  K ANSAC  I  IONS  ON  SONIC'S  AND  L’L  I  K ASONICS.  VOL  SU  26.  NO.  4.  IULY  1974 


JO? 


there  exist  one  real,  and  two  complex  conjugate  roots.  For  all 
real  media  the  Poisson  ratio  is  bounded  by 

0  <  o  <  0.5.  (S) 

Complex  roots  are  not  acceptable  for  they  imply  attenua¬ 
tion  or  damping  which  is  not  the  case  for  the  free  surface  here 
Further  the  root  V  must  have  a  magnitude  less  than  the  veloc¬ 
ity  of  the  bulk  shear  wave  in  the  medium, 

V<VS.  (6) 

The  constraint  is  a  manifestation  of  the  radiation  condition. 
Thus  there  is  only  one  root  [2)  satisfying  (6).  The  root,  given 
by 

V=VR  '  (7) 

corresponds  to  the  Rayleigh  surface  wave  propagating  with  ve¬ 
locity  VR.  The  root  (7)  always  exists.  The  Rayleigh  velocity 
varies  approximately  from 

0.87  ys  - *  0.96  Vs  (8) 

as  the  Poisson  ratio  varies  from  0  to  0  5. 


B.  The  Leaky  Rayleigh  Wave 

Consider  the  elastic  isotropic  solid  above,  now  loaded  by  a 
liquid  half-space.  Again  one  seeks  a  solution  to  the  equations 
of  motion  in  which  the  displacement  amplitude  attenuates  ex¬ 
ponentially  normal  to  the  interface.  At  the  interface,  normal 
stresses  and  displacements  must  be  continuous.  Further  tan¬ 
gential  stresses  in  the  solid  must  vanish.  Upon  application  of 
the  boundary  conditions,  a  new  secular  equation  is  derived. 
The  secular  equation  is  given  by 

4(KS/F)2[1  -  (Fs/K)2|l'2I(Ks/KD)2  -  (Fs/K)2)1'2 
+  11  -  2(Fs/F)2IJ  +  (pt/p5)((Ks/K0)2 

-  (Ks/K)2],/2[(Ks/Kl)2  -  (VSIV)2V'17  =0  (9) 


where  VL  is  the  phase  velocity  in  the  liquid.  The  densities  Pi 
and  ps  correspond  to  those  of  the  liquid  and  solid,  respectively. 

Two  permissible  roots  exist  for  (9);  one  of  them  is  real.  This 
root  always  exists  as  in  the  Rayleigh  wave  case.  The  surface 
wave  associated  with  this  solution  is  known  as  the  Stoneley 
wave.  It  is  a  pure  surface  wave  propagating  parallel  to  the  in¬ 
terface  without  attenuation  and  exponentially  damped  in  both 
directions  normal  to  the  boundary.  The  Stoneley  phase  veloc¬ 
ity  Vst  is  constrained  to  be 

VST<Vl  (10) 

There  is  a  second  root  which  is  complex  [3],  The  value  of  this 
complex  root  approaches  the  Rayleigh  wave  velocity  as  the 
density  of  the  liquid  approaches  zero.  The  inhomogeneous 
wave  associated  with  the  root  is  the  leaky  Rayleigh  wave. 

The  complex  root  is  an  allowed  solution  for  a  liquid/elastic 
solid  system.  Attentuation  implied  by  the  complex  nature  of 
the  root  has  physical  meaning  in  this  instance.  As  the  leaky 
Rayleigh  wave  propagates,  energy ’s  transferred  into  the  liquid. 
The  real  and  the  imaginary  parts  of  the  complex  root, 


Vrr  «  Re  V 
Vri  -  Im  V 


(11) 


yield  the  phase  velocity  of  the  leaky  Rayleigh  wave  and  the 
attenuation,  respectively. 


C.  Conditions  for  existence  oj  Leaky  Rayleigh  Wave 

In  contrast  to  the  free  Rayleigh  wave,  the  leaky  Rayleigh 
wave  does  not  exist  for  all  liquid/elastic  solid  plane  half-spaces. 
A  characteristic  of  the  leaky  Rajleigh  velocity  is 

1'k  ~  Vrr  <  t’s  <  l  u  (1  2) 

The  value  of  1'##  must  be  less  than  the  bulk  shear  velocity 
due  to  the  radiation  condition  Due  to  the  loading  by  the  liq¬ 
uid,  the  leaky  Rayleigh  wave  velocity  must  be  greater  than  the 
free  Rayleigh  wave  velocity. 

If  the  sonic  velocity  in  the  liquid  exceeds  the  bulk  shear 
wave  velocity,  then  no  leaky  Rayleigh  wave  can  be  propagated. 
From  consideration  of  the  radiation  condition  again,  the  leaky 
Rayleigh  wave  velocity  is  constrained  to 

VL  <  Vrr  _  U3) 

Thus,  from  (12)  and  (13).  one  obtains  a  necessary  condition 
for  existence  of  the  wave,  that  is, 

vL  <  ys  04) 

This  condition  follows  directly  from  Snell's  law:  for  a  sur¬ 
face  wave  radiating  into  the  liquid,  the  direction  of  the  radia¬ 
tion  is  given  by 

sin  dL  =  VLIVRR  (15) 

where  for  the  surface  wave,  sin  dRR  =  1 .  If  the  liquid  velocity 
exceeds  the  shear  wave  velocity,  this  implies  that  VL  >  VRR 
or  considering  (15), 

sin  0i>  \.  (16) 

The  necessary  condition  ( 1 4)  is  well  known ;  however,  it  is 
not  the  only  condition  necessary  in  order  to  admit  a  leaky 
Rayleigh  wave  solution.  A  second  condition  which  encom¬ 
passes  the  first  is  obtained  through  numerical  analysis  of  (9). 

III.  Numerical  Solution  to  Secular  Equation 

The  secular  equation  (9)  may  be  written  as 

4(Ki/K)M(Ki/Fs)2  -  (Ki/nM^KK^/Ko)2 

-  (*i/n2]'/J  +  I (VL/VS)2  -  2(Kl/K)2]2 

+  (VLIVsf(Pl/psH(VL/VDf 

-  (^t/^)2l,/2[l  -  (Vl/V)2]-'/2  =0  (17) 

for  ease  in  calculation  and  consideration  of  the  influence  of 
the  velocity  and  density  parameters.  Solutions  to  the  leaky 
Rayleigh  wave  equation,  (17),  are  obtained  numerically.  As  a 
result  of  the  calculation,  it  is  seen  that  the  existence  of  a  leaky 
Rayleigh  wave  is  dependent  upon  the  ratios:  VLIVD ,  VL/VS, 
and  Pi  IPs-  For  given  parameters  VlIVd  and  Vi/Vs.  the 
maximum  allowable  Pilpg  is  calculated.  The  results  are  tabu¬ 
lated  in  Table  I.  The  limits  on  the  Poisson  ratio  bound  the  re¬ 
lationship  of  the  bulk  shear  wave  velocity  to  the  bulk  elastic 
longitudinal  wave  velocity.  Since  the  Poisson  ratio  has  an 
upper  limit  of  1  / 2,  the  elastic  bulk  shear  velocity,  V§,  is  given 
by 

Vs<yDN 7.  (18) 

Due  to  this  constraint  the  range  of  is  limited  by  the 

parameter  VLIVD.  This  is  expressed  in  Table  I  by  the  cutoff 
in  calculated  maxima  in  the  density  ratio. 

As  can  be  seen  from  Table  I,  as  the  bulk  longitudinal  wave 
velocity  increases  for  a  given  shear  wave  and  liquid  velocity, 
the  maximum  density  ratio  increases.  Likewise,  as  the  elastic 
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I  ABLE  I 

M  Al  InUMII  I  ,)/  P\  Xs  X  Kl  M  I U  IN  III  I  /  I/, 

(Hoki/iimxi  Hi  xiiim.)  xm>  Ij  FsIVikikxi  Hi  xiiim.) 


0.05 

0.10 

0.15 

0.20 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.10 

7.56 

0.15 

4.87 

5.36 

0.20 

3.60 

3.74 

0.25 

2.83 

2.89 

3.05 

0.30 

2.30 

2.34 

2.43 

2.58 

0.35 

1.93 

1.96 

2.01 

2.09 

0.40 

1.65 

1.66 

1.70 

1.75 

1.83 

0.45 

1.42 

1.43 

1.46 

1.49 

1.55 

1.60 

0.50 

1.24 

1.25 

1.26 

1.29 

1.32 

1.37 

0.55 

1.08 

1.09 

1.10 

1.12 

1.14 

1.18 

1.21 

0.60 

0.95 

0.95 

0.96 

0.97 

0.97 

1.02 

1.05 

1.06 

0.65 

0.83 

0.83 

0.84 

0.85 

0.86 

0.88 

0.90 

0.92 

0.70 

0.72 

0.72 

0.73 

0.73 

0.74 

0.75 

0.77 

0.  79 

0.80 

0.75 

0.62 

0.62 

0.62 

0.63 

0.63 

0.64 

0 .  t>6 

0.67 

0.69 

0.70 

o.so 

0.52 

0.52 

0.53 

0.53 

0.53 

0.54 

0.55 

0.56 

0.57 

0.59 

0.85 

0.43 

0.43 

0.44 

0.44 

0.44 

0.45 

0.45 

0.46 

0.47 

0.48 

0.50 

0.90 

0.30 

0.30 

0.30 

0.31 

0.31 

0.32 

0.32 

0.33 

0.33 

0.34 

0.35 

0.95 

0.22 

0.22 

0.23 

0.23 

0.24 

0.24 

0.25 

0.26 

0.27 

0.27 

0.29 

V./v. 


TABI.F  II 

FXXVIIM  |x  III  I.IQI  III  Fl  XXlk  Sill  III  CoMHIX  Mins'.  UHII  II  F  Ml  nil 

Exixll  M  I  C'oMllllliN 


Llquld/ClastU  Solid 

V  /V 

V  D 

Vvs 

Vs 

Hg/Al 

0.226 

0.^77 

5.0 

CCl^/Lucite 

0.34b 

0.842 

1.35 

Hg/Ni 

0.265 

0.<*o5 

1.53 

Hg/Fused  Silica 

u.24i 

(  .  365 

6 . 1*. 

HjO/ Ice4 

0 .  i*29 

0.966 

1.11 

IV.  Conclusion 

The  leaky  Rayleigh  wave,  which  is  the  complex  solution  to 
(9),  approaches  the  Rayleigh  wave  as  the  density  of  the  liquid 
goes  to  zero.  While  the  free  Rayleigh  wave  always  exists,  the 
leaky  Rayleigh  wave  does  not.  It  is  well  known  that  if  the 
liquid  velocity  is  greater  than  the  bulk  elastic  shear  wave  ve¬ 
locity,  then  there  is  no  admissible  leaky  Rayleigh  solution. 
However,  as  a  result  of  parametric  investigation  of  the  secular 
equation  via  numerical  analysis  a  second  set  of  conditions  is 
evident.  Existence  of  a  solution  does  not  only  depend  on 
VLtVs  but  on  the  values  of  the  three  ratios:  Vj JVD,  VjJVs. 
and  PL/Ps- 
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Fig.  1.  Combinations  of  F^/Fy,  F iJVp.  and  p ijpy  which  allow  prop¬ 
agation  of  a  leaky  Rayleigh  wave.  Points  must  lie  within  the  volume 
shown  for  wave  to  exist. 


shear  velocity  increases  for  given  V L  and  V  D.  the  maximum 
density  ratio  increases.  Further,  for  constant  VD  and  Ky,  i.e., 
for  a  given  solid,  an  increase  in  the  liquid  velocity  decreases 
the  maximum  allowable  density  ratio. 

Fig.  1  shows  the  range  of  admissible  parameter  sets.  All  sets 
of  the  parameters  under  the  curved  surface  yield  a  valid  solu¬ 
tion  to  (1 7).  However,  if  for  a  given  set  of  Vl!Vd  and  VLIVS 
the  density  ratio  lies  above  the  curve,  then  no  leaky  Rayleigh 
wave  is  possible.  Thus  one  obtains  the  condition  for  existence 
of  the  leaky  Rayleigh  wave. 

To  illustrate  the  condition,  several  examples  are  given  in 
Table  II.  For  real  media,  there  are  cases  where  the  velocity 
and  density  ratio  parameters  are  such  that  no  leaky  Rayleigu 
propagation  is  possible.  It  is  evident  from  Table  II  that  there 
are  cases  where  a  leaky  Rayleigh  waves  cannot  be  generated 
even  if  the  usual  condition  (VL  <  Vs)  is  satisfied. 
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ACOUSTO-OPTIC  INTERACTION  OF  SECOND  HARMONICS  IN  LAMB  WAVES 
N.G.  BROWER  and  W.G.  MAYER 


Physios  Department ,  Georgetoun  University ,  Washington,  D.C.  20007 


Resume.  -  La  possibilite  de  generation  d ' harmoniques  d'ondes  de  Lambse  propageant  dans  une  plaque 
isotrope  est  examine.  Les  conditions  de  la  generation  d'harmoniques  sont  discutees.  Un  systeme  de 
detection  utilisant  1 'interaction  acousto-optique  permet  de  mettre  en  evidence  1 'existence  de  deu- 
xiemes  harmoniques  en  mesurant  les  disymdtries  de  la  figure  de  diffraction.  Les  resultats  expdrimen- 
taux  sont  presentes. 

Abstract.  -  The  possibility  of  harmonic  generation  in  Lamb  waves  propagating  on  isotropic  plates  is 
investigatea.  The  conditions  for  harmonic  generation  are  discussed.  A  detection  scheme  using  the  acous-' 
to-optic  interaction  is  designed  by  which  asymmetries  in  the  diffraction  pattern  are  measured  to  de¬ 
termine  the  existence  of  second  harmonics.  The  results  of  the  experiments  are  presented. 


1.  Introduction.  -  Finite  amplitude  mechanical  vi¬ 
brations  in  an  elastic  solid  are  nonlinear  in  natu¬ 
re  which  may  give  rise  to  a  number  of  interesting 
interactions.  One  such  interaction  is  the  generation 
of  harmonics  as  a  finite  amplitude  acoustic  wave  is 
propagating  in  an  elastic  solid. 

Acoustic  bulk  second  harmonic  generation  was 
observed  by  Gedroits  and  Krasil  'nikov^When  an  ini¬ 
tially  monochromatic  longitudinal  ultrasonic  bulk 
wave  was  launched  In  an  elastic  medium,  a  secondary 
longitudinal  wave  was  observed.  The  second  harmonic 
amplitude  increased  linearly  with  the  interaction 
length,  i.e.,  with  the  distance  of  travel. 

Harmonic  generation  of  acoustic  surface  waves 
in  crystals  was  observed  by  LPpen  /2/.  His  experi¬ 
mental  results  agreed  with  the  theoretical  analysis. 
As  in  the  case  of  bulk  waves,  the  amplitude  of  the 
second  harmonic  surface  wave  varied  linearly  with 
the  Interaction  length  and  quadratically  with  the 
fundamental  surface  wave  amplitude. 

The  possibility  of  harmonic  generation  In  elas¬ 
tic  isotropic  plates  is  Investigated  In  the  present 
experiment.  As  an  Introduction,  general  features 
of  second  harmonic  generation  are  presented  as  a 
theoretical  background.  Next,  conditions  for  gene¬ 
ration  are  discussed.  The  detection  scheme  utili¬ 
zing  a  laser  optical  probe  Is  outlined.  A  brief  des¬ 
cription  of  the  experimental  design  Is  given.  Final¬ 
ly  results  of  the  experiment  are  presented. 

2.  Theoretical  Background.  -  Using  a  perturbative 
approximation  /3/,  the  nonlinear  equations  of  mo¬ 
tion  for  an  isotropic  elastic  solid  bulk  medium  re¬ 


duce  to 

-  cf  ^  (f  .  ij)  +  C|  ?  x  ?  x  U  =  Q,  (1) 

L 

where  U  is  the  particle  displacement  amplitude  and 
CL  and  C-j.  are  the  longitudinal  and  transverse  bulk 
wave  velocities.  The  term  on  the  righthand  side  of 
eq.(l),  Q,  is  the  source  or  forcing  term.  In  the 
case  of  harmonic  generation,  Q  consists  of  terms 
quadratic  in  the  fundamental  plane  wave  amplitude. 
An  Iterative  solution/4/  to  the  nonlinear  equation 
of  motion  yields 

U2  -  C  \i\  x.  (2) 

The  amplitude  of  the  second  harmonic,  U2,  is  pro¬ 
portional  to  the  square  of  the  fundamental  wave 
amplitude,  Uj,  and  also  to  the  interaction  length, 
x,  with  C  a  proportionality  constant.  This  result 
Is  also  obtained  by  LPpen  /2/  for  surface  waves  In 
which  case  Uj  is  the  fundamental  surface  wave  am¬ 
plitude.  Similarly,  for  plate  mode  Interactions,  eq. 
(2)  Is  valid,  where  Uj  Is  the  Lamb  mode  fundamental 
amplitude.  The  characteristics  Implied  by  eq.  (2) 
will  be  used  to  Identify  a  generated  harmonic. 

In  order  that  a  harmonic  be  generated,  the  ve¬ 
locities  of  the  fundamental  and  the  harmonic  must 
be  equal,  that  Is 


C8-176 


JOURNAL  DE  PHYSIQUE 


vtuj)  =  v(2u,j).  (3) 


This  is  a  resonance  condition.  Two  means  of  satis¬ 
fying  eq.  (3)  are  possible  for  Lamb  mode  interac¬ 
tion.  The  curves  in  Fig.  1  show  the  Lamb  mode  velo- 


Fig.  1.  -  Velocity  disnersioo  curves  for  a  free 
brass  plate  as  a  function  of  product  frequency  times 
plate  thickness  (fd  in  units  of  Wiz.mm).  Selected 
points  for  harmonic  generation  are  indicated  by 
numbers  1  to  4  -  see  text. 

city  dispersion  for  an  unloaded  brass  plate.  One 
means  of  satisfying  eq.  (3)  is  by  exciting  a  locally 
non-dispersive  Lamb  mode.  For  example,  if  the  funda¬ 
mental  were  excited  at  7  Miz  in  the  Aj  mode  (point  1) 
then  a  second  harmonic  at  14  MHz  (point  Z)  will  sa¬ 
tisfy  the  resonance  condition.  The  resonance  condi¬ 
tion  may  also  be  satisfied  for  dispersive  regions 
provided  there  is  mode  coupling.  For  example,  con¬ 
sider  a  fundamental  in  the  Aj  mode  at  5  MHz  (point  3) 
the  second  harmonic  is  in  the  Aj  mode  (point  4). 
Whenever  Lamb  modes  couple,  eq.  (3)  is  satisfied. 

This  coupling  has  been  observed  by  Brower  and 
Mayer  /5/  for  the  "Isotropic  three-phonon  interac¬ 
tion". 

The  two  cases  of  possible  resonance  are  experi¬ 
mentally  investigated  here.  Table  1  lists  the  funda¬ 
mentals  used  and  the  possible  harmonics  generated. 


3.  Optical  Detection  of  Second  Harmonics.  - 
A.  Basic_Iheory . 

The  optical  detection  system  is  based  on  the 
fact  that  a  light  beam  (laser)  will  be  diffracted 
upon  reflection  from  a  periodic  surface  corugation. 
This  corrugation  may  be  an  acoustic  surface  wave 
(a  Rayleigh  wave)  or  a  plate  vibrational  mode  (a 
Lamb  wave).  A  theoretical  treatment  of  this  acousto¬ 
optic  interaction  for  a  Rayleigh  wave  has  been  gi¬ 
ven  by  Mayer,  Lamers,  and  Auth  /6 /.  This  treatment 
can  be  extended  to  include  Lamb  waves  / 7/ .  Fig.  2 
is  a  schematic  diagram  of  the  laser-surface  wave 


Fig.  2.  -  Schematic  diagram  for  acousto-optic  in¬ 
teraction  caused  by  Lamb  waves.  K  refers  to  the 
wave  vector  of  the  light  incident  at  an  angle  0  , 

9  is  the  angular  spacing  of  the  first  light  dif¬ 
fraction  order,  and  x  denotes  the  direction  of 
travel  of  the  Lamb  wave  on  the  plate. 

interaction.  The  light  is  Incident  into  the  surface 
at  an  angle  0  and  is  diffracted  into  discrete  dif¬ 
fraction  orders  with  angular  spacing  6,  centered 
around  the  reflection  angle  0.  The  light  Intensity 
in  the  mth  diffraction  order  due  to  a  single- 
frequency  Lamb  mode  of  amplitude  Uj  is  given  by 

i  ±1  -  (n>*  w 

where 


TABLE  1.  Fundamental  Lamb  Modes  and  Harmonics 


Fundamental  Mode 

Second  Hanaonlc 

Aj  at  7  MHz 

A}  at  14  PWz 

Aj  at  S  Nlz 

Aj  at  10  Miz 

Yj  *  2KUj  COS  0, 

with  K  the  wavenumber  of  the  light  and  JN  the  mth 
order  Bessel  function.  It  is  evident  from  eq.  (4) 
that  the  diffraction  pattern  is  symmetric  with  res¬ 
pect  to  the  central  order  reflected  at  angle  0. 

Now  consider  a  surface  corrugation  caused  by  a 
Lamb  wavo  which.  In  addition  to  Its  fundamental 


J 


SL 


...  .{JUT 


'  -.laWB 

■MMWiBfcaaMi* 


N.6.  BROWER,  and  al. 


C8-17- 


frequency,  contains  a  second  harmonic  component  of 
amplitude  U2-  In  this  case  the  particle  displacement 
will  no  longer  be  pure  sinusoidal  and  a  distorted 
Lamb  wave  will  interact  with  the  incident  light  beam. 
The  resulting  diffraction  pattern  will  be  asymmetric 
with  respect  to  the  central  order.  The  amount  of 
asymmetry  will  depend  on  the  magnitude  of  U2  relati¬ 
ve  to  Uj. 

This  process  is  similar  to  that  described  by 
Neighbors  and  Mayer  /8/  for  nonsinusoidal  Rayleigh 
waves.  Using  the  same  approach  for  nonsinusoidal 
Lamb  waves  and  assuming  that  only  the  fundamental 
and  a  second  harmonic  component  are  present  in  the 
wave,  the  light  intensities  in  the  positive  and  ne¬ 
gative  first  orders  are  given  by 

l+r  [••■  -  J3(t1)J2(y2)-J1(y1)J1(t2)+J1{y1)J0(t2) 

-  J3(y1)Ji(y2)  ...]  2  .  (5) 

l-r  [•••  '  J3(y1)J1(y2)-o1(t1)J0(y2)-J1(y1)J1(y2) 

-  J3(yi)02(y2)  ...)  2  ,  (6) 

where  y2  =  2KU2  cos  0. 

The  series  of  products  of  Bessel  functions  in 
above  equations  can  safely  be  truncated  to  the  terms 
shown  since  both  Yj  and  y2  are  usually  small  so 
that  Insignificant  errors  In  the  values  of  I+j  and 
I_j  are  Introduced  by  neglecting  terms  In  Jr0$  with 
r,s  >  3. 

Equations  (5-6)  are  the  bases  In  the  optical 
detection  of  the  presence  of  second  harmonics.  If 
U2  *  0,  one  finds  that  *-l  >»♦!■  giving  rise  to 
an  asymmetric  diffraction  pattern.  It  should  be  no¬ 
ted  that  eqs.(5-6)  reduce  to  eq.  (4)  if  U2  «  0.  In 
this  case  all  terms  vanish  except  Jj(yj)Jq(y2)  which 
equals  J1(y1). 

This  analysis  provides  a  method  of  determining 
the  presence  or  absence  of  second  harmonics  In  a  vi¬ 
brating  solid  plate.  Probing  the  surface  with  a  light 
beam. and  measuring  the  Intensity  distribution  In  the 
diffraction  pattern  yields  Information  about  the  ma¬ 
gnitude  of  Uj  compared  to  Uj.  Moreover,  probing  the 
turf ace  at  different  distances  from  the  source  of 


the  Lamb  wave  generation,  one  can  measure  the  chan¬ 
ge  in  harmonic  content  as  a  function  of  distance. 

Due  to  the  method  of  generating  the  funda¬ 
mental  Lamb  mode,  i.e.,  by  means  of  the  "liquid 
wedge"  method  /9/,  the  possibility  exists  that 
second  harmonic  components  are  generated  in  the  li¬ 
quid  prior  to  incidence  on  the  plate  surface.  Ho¬ 
wever,  if  such  components  should  be  present  and 
should  continue  to  propagate  along  the  plate  they 
would  at  best  decrease  in  amplitude  as  a  function 
of  distance  traveled  on  the  plate.  Their  amplitude 
could  not  increase  as  a  function  of  distance  on 
the  plate  unless  harmonic  generation  on  the  plate 
is  possible.  Thus  eq.  (2)  will  be  satisfied  only 
if  harmonic  generation  on  the  plate  itself  occurs, 
if  harmonics  propagate  on  the  plate,  asymmetry  in 
the  diffraction  pattern  will  be  observed;  to  ob¬ 
tain  numerical  values,  a  new  parameter,  lp,  is  in¬ 
troduced,  defined  as 

!-l  ■  !+l 

Ip  S  — - - ii  x  100  (7) 

!-l 

Using  eq.  (2),  it  may  be  shown  that  Ip  is  a 
linear  function  of  distance  from  the  source  of  the 
Lamb  wave  /7/  if  the  plate  is  considered  to  be  a 
dissipationless  medium.  Therefore,  Ip  can  be  measu¬ 
red  as  a  function  of  distance. 


Fig.  3.  -  Experimental  setup  for  optic  probing  of 
haneonlc  generation  In  Lamb  waves. 
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B.  Experimental  Design  and  Results. 

The  design  for  the  optical  probing  system  is 
shown  in  Fig. 3.  Light  from  a  He-Ne  laser  is  expan¬ 
ded  and  collimated  (details  of  the  necessary  optical 
components  not  shown  in  Fig.  3.)  The  collimated 
light  beam  is  reflected  from  a  highly  polished  brass 
plate  which  supports  the  Lamb  wave.  The  diffracted- 
reflected  light  is  focussed  approximately  8  m  from 
the  plate.  This  distance  is  sufficient  to  properly 
separate  the  diffraction  orders  so  that  intensity 
measurements  of  I+j  and  I_^  can  be  made  with  a  pho¬ 
todiode. 

Measurements  of  Ip  were  made  at  various  dis¬ 
tances  x  on  the  plate  for  the  different  fundamental 
frequencies  and  modes  listed  in  Table  I.  Excitation 
of  the  modes  was  accomplished  by  adjusting  the 
quartz  transducer  to  the  required  angular  position 
for  the  generation  of  the  desired  Lamb  mode. 

The  polished  brass  plate,  excited  in  the 
mode  at  7  MHz,  was  probed  with  the  He-Ne  laser.  A 
diffraction  pattern  was  observed  and  the  asymmetry 
between  the  first  diffraction  orders  was  measured.  At 
a  relative  interaction  length  of  0.8  cm,  the  percent 
asymmetry  was  approximately  15%.  The  laser  probe  was 
then  moved  along  the  plate  in  the  positive  x-direc- 
tion  and  the  asymmetry  increased.  At  a  distance  of 
1.2  cm,  a  29%  asymmetry  was  observed.  Likewise,  for 
the  plate  in  the  fundamental  A^  mode,  excited  at 
5  MHz,  asymnetry  was  observed.  In  this  case  the 
asymmetry  also  Increased  with  interaction  distance. 
The  results  are  listed  in  Table  II. 

TABLE  II.  Ip  as  a  function  of  x. 


X 

I 

P 

0.8  cm 

15% 

1.2  cm 

29% 

Aj  mode  at  7  ttiz 

0.9  cm 

14% 

Aj  mode  at  5  MHz 

1.3  cm 

39% 

As  a  control,  the  $q  mode  at  3  MHz  was  excited 
(point  S  in  Fig.  1.).  As  can  be  seen  in  Fig.  1., 
there  Is  no  second  harmonic  which  satisfies  Eq.  (3) 
because  the  required  location  (point  t)  does  not 
fall  on  a  curve.  Therefore,  there  should  be  no  se¬ 


cond  harmonic  content  in  the  detected  mode.  The 
measured  diffraction  pattern  showed  this  to  be  true: 
no  detectable  asymmetry  in  the  first  orders  could 
be  found. 

In  conclusion,  second  harmonic  generation  in 
Lamb  modes  was  detected.  Harmonics  were  observed 
for  both  a  locally  dispersionless  mode  and  also  for 
the  mode  coupling  case,  provided  the  appropriate 
conditions  for  harmonic  generation  in  plates  were 
satisfied. 
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